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ABSTRACT 

The Pyxis globular cluster is a recently discovered globular cluster that lies 
in the outer halo (R gc ~ 40 kpc) of the Milky Way. Pyxis lies along one of the 
proposed orbital planes of the Large Magellanic Cloud (LMC), and it has been 
proposed to be a detached LMC globular cluster captured by the Milky Way. 
We present the first measurement of the radial velocity of the Pyxis globular 
cluster based on spectra of six Pyxis giant stars. The mean heliocentric radial 
velocity is ~36 km/sec, and the corresponding velocity of Pyxis with respect to 
a stationary observer at the position of the Sun is ~ —191 km/sec. This radial 
velocity is a large enough fraction of the cluster's expected total space velocity, 
assuming that it is bound to the Milky Way, that it allows strict limits to be 
placed on the range of permissible transverse velocities that Pyxis could have 
in the case that it still shares or nearly shares an orbital pole with the LMC. 
We can rule out that Pyxis is on a near circular orbit if it is Magellanic debris, 
but we cannot rule out an eccentric orbit associated with the LMC. We have 
calculated the range of allowed proper motions for the Pyxis globular cluster 
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that result in the cluster having an orbital pole within 15° of the present orbital 
pole of the LMC and that are consistent with our measured radial velocity, 
but verification of the tidal capture hypothesis must await proper motion 
measurement from the Space Interferometry Mission or HST. A spectroscopic 
metallicity estimate of [Fe/H] = —1.4 ± 0.1 is determined for Pyxis from several 
spectra of its brightest giant; this is consistent with photometric determinations 
of the cluster metallicity from isochrone fitting. 



Subject headings: globular clusters: individual (Pyxis) 



1. Introduction 

Evidence continues to accumulate that the outermost Milky Way globular clusters may 
not have originated in the same process that formed the inner globular clusters. Based 
on the recognition that the second parameter effect of horizontal branch morphology in 
globular clusters is found predominantly among outer halo (R gc > 8 kpc) clusters, |Searle fc| 
Zinn (1978)| proposed that the outermost globular clusters may have formed in chemically 



distinct "fragments" that later fell into the Milky Way halo. Building on the suggestion by 
Kunkel fc Demers (1976)] that several red horizontal branch (second parameter) globular 



clusters were potentially associated with the Magellanic Plane group of dwarf galaxies, 
Majewski (1994)| showed that there is a planar alignment between a particular sample of 



second parameter globular clusters and the Milky Way dwarfs. Recently, |Palma et al. 



2000) reaffirmed that there may be a dynamical relationship between the second parameter 



globular clusters and the Milky Way dwarf satellites. 



The Pyxis globular cluster ( [Irwin et al. 1995| ; pa Costa 1995| ) at R gc = 41 kpc 



(Sarajedini & Geisler 1996) defines the inner edge of the prominent gap in the globular 
cluster radial distribution between 40 <^ R gc ^ 60 kpc. The presence of this gap has been 
used to argue that the primordial Galactic globular cluster system ends at ~40 kpc while 



the distant, R gc > 60 kpc clusters originated in Galactic satellite dwarf galaxies (e.g., |Zinn 
1985|) . Although Pyxis lies among the "inner group" of globular clusters (i.e., inside the 
gap in R gc ), Irwin et al. (1995) propose that Pyxis may be a captured LMC globular cluster 
based on the young age they infer for the globular cluster and on its proximity to the plane 



of the LMC orbit derived from the [Jones et al. (1994)] proper motion. Further support for 



the tidal capture hypothesis comes from [Palma et al. (2000)| , where a statistical analysis of 



the likely orbital poles of the Galactic satellite galaxies and the globular clusters identifies 
Pyxis, NGC 6229, NGC 7006, and Pal 4 as the globular clusters most likely to share a 
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common orbital pole with either the Magellanic Plane galaxies (the LMC, the SMC, Draco, 
and Ursa Minor) or the Fornax-Leo-Sculptor Stream galaxies. However, these postulations 
on the origin of Pyxis have been made without the benefit of any kinematical data on the 
cluster. 

Although deep photometry of Pyxis exists (Sarajedini & Geisler 1996), no spectroscopic 
observations have been published. Indeed, Pyxis is one of the last few known globular 
clusters lacking a radial velocity (cf. Harris 1996). We report here on du Pont 2.5-m 
Telescope spectroscopic observations of Pyxis stars (§2). With our derived radial velocity 
for the cluster, we re-address the stripped LMC hypothesis for Pyxis' origin (§3), but point 
out that, in the end, we can only make predictions on the proper motions expected under 
this scenario. Unfortunately, the proper motion is required for a definitive solution to the 
question of the cluster's origin. 

2. Observations 

On the nights of 17 January and 20 January 1997, the 2.5-m du Pont Telescope 
at the Las Campanas Observatory was used to obtain spectra of Pyxis giant stars with 
magnitudes of R ~ 18. A finding chart for these stars made from the Digitized Sky Survey 
( [Lasker et al. 1990|) is presented as Figure 1. Observations of the Pyxis stars made use 
of the modular spectrograph and the SITe2 detector with a 1200 lines/mm grating. The 
selected wavelength region with this setup was approximately 7700 - 8750 A at ~ 1.3 A per 
pixel resolution. Typical exposure times were 900 seconds per observation, which provided 
enough signal-to-noise to measure radial velocities. 

Radial velocities were measured for six Pyxis giant stars by cross-correlating their 
spectra against those of bright radial velocity standard stars. The cross-correlation peak 
yielded the adopted radial velocity difference between the Pyxis stars and the standard. 
The star HD80170, a K5 giant, was observed multiple times on both nights to serve as 
the main standard. Two other HD stars and 11 LMC globular cluster stars were also 
observed as radial velocity calibration objects. All standards with known radial velocities 
were used to remove any nightly zero point offsets in the radial velocities determined by 
cross-correlating against the HD80170 template. The Pyxis stars were observed multiple 
times, and the radial velocity for each star was determined by taking the median of the 
values for a particular star. In Table 1 we present the positions, magnitudes, and velocities 
of the stars observed. We find a mean radial velocity for the six stars of 34.3 km/sec with 
a = 4.6 km/sec. If we ignore the one outlier, Pyxis D, the mean radial velocity for the five 
remaining stars increases slightly, to 35.9 km/sec with a = 2.5 km/sec. 
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2.1. Call Triplet Metallicity Estimate 

The wavelength range of our spectra includes the Ca II triplet lines at AA 8498, 8542, 
and 8662 A. [Rutledge et al. (1997a)| measured Ca II equivalent widths for a large sample of 



globular cluster stars and in a companion paper (Rutledge et al. 1997b) derived a conversion 
from Ca II reduced equivalent width, W, to metallicity. To determine a spectroscopic 
metallicity for Pyxis, we attempted to measure the strengths of these lines in our spectra. 
Of the six Pyxis stars observed, only the brightest, Pyxis A, produced equivalent width 
measurements with errors small enough to give a reliable estimate of the metallicity. The 
technique used to determine the equivalent widths for the Pyxis A Ca II lines was nearly 
identical to that used in [llutlcdge et al. (1997a) . 



For each observation of Pyxis A, the continuum was determined by linearly interpolating 
the average intensity in the |Rutledge et al. (1997a)| continuum bandpasses. The equivalent 



width was then calculated by integrating the difference between the fit continuum and 
the line feature over the line bandpass. The line feature was fit with a Gaussian function, 
and the integral was performed using the fit rather than numerically integrating the data 
since [Rutledge et al. (1997a)| concluded that this technique is preferable to direct numerical 
integration. The three individual lines were combined into a single index, ECa, following 
the method of |Rutledge et al. (1997a)|, where ECa = 0.5A 84 9 8 + A 85 4 2 + 0.6A 8662 . For Pyxis 



A, the result is ECa = 4.5 ± 0.2 A. 



The Rutledge et al. (1997a)| method for converting ECa to metallicity requires the 



reduction of the equivalent width to the value for giants at the level of the horizontal branch, 
W, so that a mean value for all cluster stars can be obtained. This is done by adopting a 
slope A(ECa)/A(VH\B — V) and extrapolating the calculated width to the expected value 
at the magnitude of the horizontal branch, Vhb- All of the published photometry for Pyxis 
stars has used the B and R bands, so we have had to approximate A(Vhb — V) using 
color-color relations for giant stars. Adopting the photometry from [Sarajedini &; Gcislci 



1996), we find that Pyxis A has R = 17.08 and that the horizontal branch of Pyxis is at 



R = 18.75. Caldwell et al. ( 1993 )| find an almost linear relationship between (B — V) and 



(V — R) for giant stars, so we have determined rough V magnitudes for Pyxis A and the 
Pyxis horizontal branch stars by estimating a (V — R) color from the (B — R) color given in 
|Sarajedini fc Geisler (1996)j We estimate that for Pyxis, Vhb = 19.25 and Vp yxis A = 17.77, 
or A(Vhb — V) = 1.48. In the table of globular cluster properties by [Harris ( 1996 )| , the 
magnitude of the horizontal branch is also given as Vhb = 19.25, so the adopted (V — R) 
colors are most likely a good approximation to the true colors. 

The reduced Ca II equivalent width of Pyxis A is therefore W = 3.6 if we follow 
Rutledge et al. ( 1997a )| and adopt a slope of A(ECa)/A(V#_B — V) = 0.62 A/magnitude. 



- 5 - 



This value is simply an estimate, since there is a dispersion of 0.2 A in the values of ECa 
from the four individual observations and since there is some uncertainty in A(Vhb ~ V)i 
probably of order 0.1 magnitudes. However, we can use this determination to get an 
estimate of the metallicity of Pyxis A for comparison with the photometrically determined 
metallicity estimates of Irwin et al. (1995) and |Sarajedini fc Geisler (1996)| . Using the 
IRutledge et al. ( 1997b )| calibrat ion of W to Zinn-West metallicity, the reduced equivalent 
width measured for Pyxis A implies a metallicity of [Fe/H]^vK = —1.4 ± 0.1. This value 
for Pyxis A is more metal-poor than the photometrically derived values of |Sarajedini fc| 
Geisler (1996)| and Irwin et al. (1995), who estimate — 1.2 =h 0.15 and —1.1 ±0.3 respectively, 



however it is consistent within the overlap of the la error bars. Any systematic error that 
leads to an underestimated equivalent width for Pyxis A results in a smaller determined 
metallicity. An error of 10% in the equivalent width measured for Pyxis A is enough to 
bring the metallicity up to —1.2 and into better agreement with the photometric values. If 
the equivalent width measurement is correct, then it is unlikely that the metallicity is much 
higher than [Fe/Hj^vy = — 1-4, since an unlikely error in the V magnitude of Pyxis A of 0.7 
magnitudes is required to raise the metallicity of the star to —1.2. 



3. Discussion 

The observations presented here were partially motivated by the possibility that the 
Pyxis globular cluster was captured from the LMC by the Milky Way. This assertion was 
originally made by Irwin et al. (1995), who noted that Pyxis, at (/,&) = (261.3,7.0)°, lies 
within a few degrees of the orbital plane of the LMC determined from the |Jones et al. (1994) 



proper motion. Further support for this hypothesis is provided by Palma et al. (200~0]| , who 



place Pyxis in a group with Pal 4, NGC 6229, and NGC 7006 as the most likely globular 
clusters to share a common orbital pole with the Magellanic Plane galaxies (the LMC, 
the SMC, Ursa Minor, and Draco). Although full space motion information is required to 
verify the Irwin et al. (1995) hypothesis, a radial velocity can provide some constraints on 
the shapes of allowed orbits for a cluster if the magnitude of the radial component is a 
significant fraction of the expected magnitude of the space velocity. 

The essence of the argument given in Palma et al. (2000) to support a capture origin 
for the Pyxis globular cluster is as follows: if it is assumed that Pyxis was captured recently 
from the LMC by the Milky Way, then the orbital pole of Pyxis is likely to be aligned with 
that of the LMC (both the LMC and Pyxis are far enough from the Galactic Center that 
precession will not significantly affect the positions of their orbital poles over a Hubble 
time). If one assumes rotation about the Galactic Center, the direction of the orbital pole of 
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the LMC can be determined by taking the cross product of the Galactocentric radius vector 
to the LMC and its space motion vector. If one accounts for the space velocity vector error, 
the position of the orbital pole of the LMC can only be confined to a family of poles along 
an arc segment in Galactocentric coordinates (cf. Figure 1 in Palma et al. 2000). Since 
the space motion of Pyxis is currently unknown, its orbital pole is not well constrained. 
However, the orbital pole can be assumed to be perpendicular to its current Galactocentric 
position, so the direction of Pyxis' orbital pole should lie on the great circle that contains all 
possible normals to its current radius vector (cf. [Lynden-Bell fc Lynden-Bell 1995|) . Figure 



2 shows an Aitoff projection of the sky in Galactocentric coordinates. The arc segment that 
defines the possible locations of the LMC's orbital pole (based on the Jones et al. [1994] 
proper motion, as adopted by Palma et al. [2000]) is shown as well as the great circle along 
which lies all possible orbital poles of Pyxis. That these two families of possible orbital 
poles for the LMC and for Pyxis intersect (at Galactocentric (l,b) = (163, —22)°) indicates 
that it is possible for these two objects to share a coplanar orbit with a common direction 
of angular momentum. It now remains to be seen whether our derived radial velocity can 
clarify whether this is likely, i.e., is the orbital pole of Pyxis likely to be near the crossing 
point of the two orbital pole families? 

We start our analysis by adopting a simple strawman model wherein Pyxis is following 
a circular orbit that is nearly coplanar with the orbit of the LMC. Simulations of tidal 
stripping of dwarf galaxies by the Milky Way (Johnston 1998) show that the debris is 
distributed around the orbit of the parent satellite with a spread in energy given by 

d$ / m sat \ 1/3 2 n 2 /x 
AE = T ^dR W UcJ *W = AW (1) 



where r ti d e is the tidal radius of the satellite, $ is the parent galaxy gravitational potential, 
f circ is the circular velocity of the Galactic halo, m sat is the satellite's mass, M Ga i is the 
mass of the parent galaxy enclosed within the satellite's orbit, and the last equality defines 
the tidal scale f. Thus, the spread in energy translates into a characteristic angular width / 
(in radians) to the debris. Taking reasonable values for the mass of the LMC and the Milky 
Way, the value of / for any debris pulled from the LMC corresponds to roughly 15°. In 
Figure 2, an arc segment along the great circle of possible orbital poles for Pyxis is marked; 
this arc segment is defined by a length within ±15 degrees of the intersection point with the 
possible orbital pole of the LMC, and indicates expectations for the orbital poles of LMC 
debris at the position of Pyxis on the sky. 

It is straightforward to derive the direction of the space motion vector required for 
Pyxis to follow a circular orbit and have an orbital pole along the arc segment in Figure 
2, i.e. in the direction of P, or (I, b) = (163, —22)°. The geometry is illustrated in Figure 
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3. The Galactocentric, Cartesian radius vector of Pyxis is (X,Y,Z) = (—13.9, —38.6,4.8) 
kpc (assuming that the distance from the Sun to the Galactic Center, R = 8 kpc, and 
where Xq = —8.0 kpc). The unit vector P from the Galactic center in the direction of 
the orbital pole at Galactocentric (/, b) = (163,-22)° is (X, Y, Z) = (-0.89,0.27,-0.37). 
In Figure 3, this vector has been translated to the location of Pyxis. The vector that is 
mutually perpendicular to the Galactocentric radius vector of Pyxis and to the orbital 
pole P gives the direction of the space motion of the Pyxis globular cluster for a circular 
orbit around the Galactic center. The unit vector direction of this space motion is 
Vdrc — {Vx, Vy, Vz) = (—0.32, 0.23, 0.92). Since our line of sight to Pyxis is mostly in the 
— Y direction and since V C i rc is mostly in the +Z direction, clearly any component of V C i rc 
along our line of sight will be small. Adopting the basic solar motion of (9, 12, 6) km/sec 
( [Mihalas fc Binney 1981| ) and a rotational velocity of the local standard of rest (LSR) of 



220 km/sec, then the component of the Sun's velocity along the line of sight from the Sun 
to Pyxis is vq los ~ —227 km/sec (the negative sign here indicates that the component of 
the Sun's velocity with respect to the Galactic Standard of Rest along the line of sight to 
Pyxis is in the sense of receding from Pyxis; see the Appendix for a discussion of the sign 
conventions used in reducing the radial velocity to a vgsr)- Since the heliocentric radial 
velocity measured for Pyxis is the difference between the intrinsic radial velocity of Pyxis 
with respect to the position of the Sun, vqsr, and the magnitude of the Sun's velocity 
projected along the line of sight to Pyxis {vheiio = vgsr, ~ v Qlos)i the globular cluster 
would have a large, positive heliocentric radial velocity if it were following a circular orbit 
in the plane defined by the pole at P. 

It is the solar motion that dominates the radial velocity in the circular orbit 
case. For example, if we assume that Pyxis has a velocity that is approximately 
the circular velocity of the Galaxy at 40 kpc, or ~ 200 km/sec, its space 
motion would then be (Vx,Vy,Vz) = (—64,46,184) km/sec. Since the Galactic 
Cartesian unit vector in the direction of Pyxis from the Sun (the line of sight) is 
(X,Y,Z) = (—0.15,-0.98,0.12), the component of this space velocity along the line of 
sightf] would be (-64 x -0.15) + (46 x -0.98) + (184 x 0.12) or v GSR = -33 km/sec (the 
direction and magnitude of this component of the vqsr i n the circular orbit case is shown 
as Vc in Figure 3, while the direction and magnitude of the measured vgsr for Pyxis is 
shown as Vm)- Including the component of the Sun's velocity along the line of sight (—227 
km/sec), the heliocentric radial velocity (Vcheiio in Figure 3) for Pyxis in this case would be 
—33 — (—227) = 194 km/sec. Since this is much larger than the measured value (36 km/sec, 



3 Again the negative sign is in the radial system of the Sun, and it indicates that this velocity points 
towards the Sun. In the Galactic Cartesian system this velocity is positive in the X, Y, and Z directions. 
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VMheiio in Figure 3), we conclude that Pyxis is not in an orbit that gives rise to a present 
space velocity near V C i rc and, therefore, the strawman model of Pyxis being on a circular 
orbit and sharing the LMC orbital pole grossly fails expectations. Thus, one or both of the 
assumptions in the strawman model must be invalid: either Pyxis is not on a circular orbit 
or/and Pyxis' orbit does not share a pole with the LMC. 

Since the inferred radial velocity of the Pyxis globular cluster with respect to a 
stationary observer at the location of the Sun has a large magnitude, Vqsr ~ —191 km/sec, 
constraints can be placed on non-circular orbits Pyxis may follow that also share the plane 
and direction of rotation of the LMC's orbit. We have calculated the orbital pole for the 
Pyxis globular cluster given all possible, realistic proper motions and accounting for our 
derived radial velocity. In Figure 4, we present the region in proper motion space that 
produces an orbital pole for the Pyxis globular cluster that is within 15° of the pole of the 
LMC's orbit. We limit the possible proper motions to those that yield a space velocity less 
than the escape velocity from the Milky Way at the position of Pyxis (~415 km/sec in the 
Galactic model of Kochanek [1996]) and obtain the shaded region in Figure 4. This is a 
prediction for the magnitude and direction of the proper motion of Pyxis with respect to 
the Sun assuming that the Milky Way capture from the Magellanic Clouds hypothesis is 
correct. 

The proper motions in the shaded region of Figure 4 are those that can produce an 
orbital pole in the direction of P given a vqsr of —191 km/sec for Pyxis. We have calculated 
the shapes and energies of the orbits allowed for Pyxis to determine if a pole in the direction 
of P is only likely for a very restricted range of conditions. For example, is the LMC capture 
origin for Pyxis only viable if Pyxis is following an extremely eccentric orbit? In fact, a 
range of orbits is possible for Pyxis given a proper motion in the shaded region in Figure 4. 
For a given elliptical orbit, the angle between the instantaneous velocity and radius vectors 
varies with position along the ellipse and that angle has a well defined minimum value for 
a given orbital eccentricity. For each space motion derived from our radial velocity and a 
proper motion from Figure 4, we have determined the angle between the velocity vector and 
the present Galactocentric radius vector for Pyxis. Assuming a closed, elliptical orbit, we 
can determine the lower limit for the eccentricity of the associated orbit having the given 
angle between the velocity and radius vectors at the present position of Pyxis. The orbits 
determined for Pyxis given our radial velocity and a proper motion in the shaded region in 
Figure 4 have eccentricities of e > 0.70, with the peak of the distribution of all allowable 
eccentricities near e ~0.8. 



We note here that few of the globular clusters with measured proper motions are 
following nearly circular orbits. Pinescu et al. ( 1999 )| has compiled all of the measured 
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proper motions for a sample of 38 Galactic globular clusters and integrated orbits for each 
cluster. Figure 5 presents a histogram of the orbital eccentricities that pinescu et al. 



(1999)| calculated for the globular clusters in their sample. The open histogram in Figure 



5 represents the data on the whole sample, while the hatched histogram represents the 
data on the 10 clusters with apoGalactica greater than 20 kpc. Less than half of the entire 
sample have eccentricities of e < 0.5, and, more importantly, for the outer halo globular 
clusters the measured eccentricities are mostly found in the range 0.6 < e < 0.8. Therefore, 
one might conclude that it is more likely than not that Pyxis is following an eccentric 
orbit and its space motion is not perpendicular to its current position. However, such a 
conclusion must be tempered with the acknowledgement of a potential selection bias for the 
latter subsample. The majority of the globular clusters with measured proper motions are 
those that are currently close to the Sun. Therefore, the clusters with large apoGalactica 
that have measured proper motions are, for the most part, currently near periGalacticon 
and thus must be following eccentric orbits. Thus, the sample of globular clusters that 
make up the hatched histogram in Figure 5 may be selected preferentially from the sample 
of outer halo globular clusters on eccentric orbits. Since few outer halo globular clusters 
near apoGalacticon have measured proper motions, the true distribution of eccentricities 
for outer halo globular clusters is unknown. However, the fact remains that a non-negligible 
(and perhaps dominant) fraction of the outer halo globular cluster population is orbiting 
the galaxy with eccentricities near 0.8, and since the majority of the outer halo globular 
clusters with orbits integrated by pinescu et al. (1999)| have eccentricities near e ~0.8, it is 



at least conceivable that the orbit of Pyxis has a similar eccentricity. 

For our measured Pyxis radial velocity, we have integrated orbits for a grid of > 1500 
proper motions found in the shaded region in Figure 4 in the potential of Johnston et 



al. ( 1995")] for 10 Gyr each. The orbital energy of Pyxis determined from the majority of 



the proper motions that produce an orbital pole at P is within the 1 — o error bars of 
the orbital energy of the LMC, although the error bar is large (Elmc — —2.1 ± 0.9 x 10 4 
km 2 /sec 2 ). However, Johnston (1998) found that in simulations of tidal stripping, debris 
was found within ±3AE of the parent object (see eq. 1). All of the orbits produced from 
our measured radial velocity and a proper motion in the shaded region of Figure 4 are 
within ±3AE of the orbital energy of the LMC. Since the P orbits do not require extremely 
unlikely constraints on the eccentricity and since the orbital energies for these orbits are 
similar to expectations for debris from the LMC, the LMC capture origin for Pyxis remains 
viable. 
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4. Conclusions 

It has been proposed since its discovery that the Pyxis globular cluster may have been 
captured by the Milky Way from the Magellanic Clouds. If the space motion for Pyxis 
were known, a comparison of the position of its orbital pole with respect to the LMC as 
well as a comparison of its angular momentum and orbital energy to that of the LMC 
would allow one to determine if the two objects share similar orbits. Although only one 
component of the space motion of Pyxis is now measured, some constraints can be placed 
on its possible orbit in the tidal capture scenario. A circular orbit with an orbital pole at 
(I, b) = (163, —22)° is completely ruled out by the measured radial velocity. However, we 
have shown here that the large radial velocity of Pyxis with respect to a stationary observer 
at the position of the Sun does not rule out the possibility that the cluster was captured 
from the LMC since a reasonable range of viable orbits with e ~0.8 exist for Pyxis that are 
also similar in energy and angular momentum to that of the LMC. No suitable first epoch 
plate material is known to exist for Pyxis, so an attempt to measure its proper motion to 
better determine the likelihood that Pyxis may be a captured LMC globular cluster will 
require precise observations with the HST or the Space Interferometry Mission (SIM). 

Although proper motions are not available for the majority of the outer halo globular 
clusters, their spatial distribution has been used to argue that they are likely to have been 
accreted into the halo (e.g., Majewski 1994, Palma et al. 2000). Recently, Dinescu et al. 
(2000) have measured a proper motion for the young globular cluster Pal 12 and they 
find that its orbit is what one would expect if it had been captured from the Sagittarius 
dwarf galaxy. An accretion origin of the outer halo, second parameter horizontal branch 
globular clusters is often invoked to explain the possible younger age of some of these 
objects (where youth is inferred either from the second parameter effect itself or from 
relative age estimates determined from the cluster CMDs). The physical mechanism that 
causes the second parameter effect in globular clusters is still unknown: Although it is now 
generally agreed that there are indeed some globular clusters with anomalously young ages, 
age differences alone may not be enough to explain the second parameter effect. Whether 
or not the physical mechanism that causes the effect is age, the possibility that conditions 
somehow favor the formation of second parameter globular clusters preferentially in Milky 
Way satellite galaxies (which later get accreted by the Galaxy) may explain the source of 
the differences between second parameter and non-second parameter globular clusters. 



The age measurement for the Pyxis globular cluster by |Sarajedini fc Ceisler (1996 



13.3 ± 1.3 Gyr, suggests that it is younger by ~3 Gyr than the oldest Milky Way globular 
clusters when measured on the same age scale. Recently, age measurements for the oldest 
LMC globular clusters have been made QOlsen et al. 1998| ) using a different technique than 
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that used for Pyxis, but their average age of 15.3 ± 1.5 Gyr places them similar in age to 
the oldest Milky Way globular clusters, when calibrated onto the same absolute age scale. 
Another study of a different sample of LMC clusters QJohnson, et al. (1999)[ ) also finds the 



oldest LMC clusters to be as old as the old Milky Way clusters. Thus, we may conclude 
that typical LMC clusters are older than Pyxis. However, at least one of the clusters in the 
Olsen et al. (1998) sample is ~ 2 Gyr younger than the others (NGC 1898), which makes it 
similar in age to Pyxis. Therefore, it is not impossible to place Pyxis in the "LMC family" 
of clusters from age arguments, though it does appear that Pyxis would be at the young 
end of the age range for old LMC clusters. 

It may be noted, however, that the current orbital pole of the Small Magellanic Cloud 
(SMC) is also very near the intersection point of the poles of the LMC and Pyxis (see 
Figure 2). Since the SMC is more fragile due to its weaker gravitational potential, perhaps a 
more attractive origin for Pyxis is from stripping of the SMC rather than the LMC. Recent 
studies of SMC globular clusters have found that the SMC clusters show a range in ages 
(e.g., [Shara et al. 1998| , Mighell ct al. 1998| ) including at least one cluster with an age similar 



to Pyxis (NGC 121). The orbital energy of the SMC has a larger magnitude and a smaller 
error bar than that of the LMC, so not all of the orbits produced from a proper motion in 
the shaded region in Figure 4 have orbital energies similar to expectations of SMC debris. 
Only the orbits having proper motions found in the inner part of the shaded region, with 
a total magnitude of the proper motion of ~0.75 mas/yr, have orbital energies consistent 
with an SMC capture origin. Since the same orbital energy and age arguments applied to 
support the LMC capture origin also apply to the SMC, we consider it a possibility that 
Pyxis may have been captured from either the LMC or the SMC. 
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A. Definition and Sign Conventions for Vqsr 

The conversion of velocities among various reference frames is treated in the literature 
and in the standard texts, such as Mihalas & Binney (1981). Because we have found some 
confusing misuse of the standard terminology in the literature, we provide this detailed 
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explanation of our sign and naming conventions. 

The radial velocity that one measures for a star is the velocity of that object with 
respect to the Earth. Often, corrections are made to this velocity to remove the motions of 
the Earth and Sun, which reduces the measured radial velocity to a velocity with respect 
to some standard of rest. For example, the measured heliocentric radial velocity (vheiio) 
is reduced to the radial velocity with respect to the Local Standard of Rest (vlsr) by 
removing the Sun's peculiar velocity with respect to the LSR: 



This velocity can be further reduced to the Galactic Standard of Rest (GSR) by removing 
the Sun's orbital velocity around the Galactic Center. So, 



Referring to this velocity as "vgsr" or a "Galactocentric" velocity apparently causes 
some confusion in the interpretation of velocity data. When reduced using the above 
two equations, the velocity referred to as vqsr is the velocity of the object as seen by a 
stationary observer at the position of the Sun. The direction of this velocity is along the 
line of sight between the object and the Sun and not along the line of sight between the 
object and the Galactic Center. The latter misinterpretation of "Galactocentric" velocity 
(as we have found in some articles in the literature) can lead to misleading or erroneous 
conclusions. 

There is an additional ambiguity in the definition of vqsr, and that is the sign 
convention. For a typical radial velocity, positive refers to a velocity that is moving away 
from the origin, and negative refers to a velocity that is approaching the origin. The origin 
for the "Galactocentric" radial velocity, or vqsr, is the Sun and not the Galactic center. 
Therefore, the sign convention for vqsr is that a positive velocity indicates that the object 
is moving away from a stationary observer at the position of the Sun and a negative velocity 
indicates that the object is moving towards a stationary observer at the position of the 
Sun. The right hand sides of the two equations above are collectively the velocity of the 
Sun along the line of sight. This sign convention introduces additional confusion because 
the sign may not agree with the sign convention for the Cartesian Galactocentric (U, V, W) 
system, and because the sign of the contribution of the Sun's motion, vq los can seem 
counterintuitive. 



vlsr = Vheiio + [9cos(6) cos(Z) + 11 cos(6) sin(/) + 6 sin(fo)] km/sec. 



(Al) 



Vgsr = v L sr + [220 cos(6) sm(l)]km/ sec. 



(A2) 



For example, the Pyxis globular cluster has been measured to have a heliocentric radial 
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velocity of ~36 km/sec. Using the above equations, v G sr — — 191 km/sec for Pyxis. The 
proper interpretation of this velocity is that Pyxis is approaching a stationary observer 
at the position of the Sun with a velocity magnitude of 191 km/sec. However, the Sun is 
located in Galactic Cartesian coordinates at (X, Y, Z) = (—8, 0, 0) and Pyxis is located at 
(X,Y,Z) = (—13.9, —38.6,4.8). Therefore, the components of this velocity, vgsr — —191 
km/sec, in Galactic Cartesian coordinates are positive in X, positive in Y, and negative in 
Z. Moreover, even though the Sun's motion is increasing the separation of Pyxis from us 
(i.e., increasing the recessional velocity) vq los is negative. 
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Fig. 1. — A finding chart for the six Pyxis stars observed. The image was created using the 
Digitized Sky Survey ( |Lasker et al. 1990| ). In this image, the center is a 2 ooo.o — 9 h 07 m 57.8 s , 
£2000.0 = — 37°13'17". North is up and East is to the left. 



-17- 



NGP 








SGP 



longitude 



Fig. 2. — An Aitoff projection of the sky in Galactocentric coordinates. The solid arc near 
the center traces the locations of all possible poles for the orbit of the LMC determined from 
the lJones et al. (1994)| proper motion. The locations of the possible orbital poles for the SMC 
(using the Irwin et al. (1996) proper motion) are shown, as well. The dashed great circle is 
the family of all possible normals to the current location of the Pyxis globular cluster. The 
orbital pole of Pyxis is assumed to lie somewhere along this great circle. The solid arc along 
the great circle contains all points within 15° of the intersection between the LMC's possible 
orbital poles and those of Pyxis. If Pyxis has been recently captured from the LMC by the 
Milky Way, its orbit should share nearly the same pole as the LMC, and it is likely to lie 
along the arc shown here. 
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Fig. 3. — see next page 
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Fig. 3. — A 3-D projection of the sky in Galactic Cartesian coordinates. The Sun is located 
at (X,Y,Z) = (-8,0,0) kpc and Pyxis lies at (-13.9,-38.6,4.8) kpc. The shaded plane 
contains the line of sight connecting the Galactic Center to Pyxis and is perpendicular to 
the vector P, which is a unit vector that points in the direction of the orbital pole located 
at {l,b) = (163, —22)°, a pole that Pyxis may share with the LMC. Any velocity vector 
for Pyxis that lies in the shaded plane will give Pyxis an orbital pole in the direction of P 
(by definition). The circular motion orbital vector v circ is a unit vector that is mutually 
perpendicular to both P and the Pyxis/Galactic Center line of sight and therefore lies in 
the shaded plane. If Pyxis is following a circular orbit with its pole at (163, —22)°, then 
its space velocity should lie along v C i rc . Assuming that the space velocity of Pyxis has a 
magnitude near the circular velocity of the Galaxy at 40 kpc, or ~200 km/sec, the predicted 
radial velocity (as seen by a stationary observer at the location of the Sun) for Pyxis is V c , 
or —33 km/sec, which corresponds to a heliocentric radial velocity (V c heiio) of 194 km/sec. 
Given that the solar motion along the line of sight to Pyxis (Vs U n) is —227 km/sec and 
that the measured heliocentric radial velocity of Pyxis (V m heiio) is 36 km/sec, the measured 
radial velocity of Pyxis with respect to a stationary observer at the location of the Sun 
(V^j), is —191 km/sec, ruling out a circular orbit for Pyxis with a pole in the direction of 
P. The unknown proper motion of Pyxis is only constrained to be perpendicular to V rn (by 
definition). It is plausible that Pyxis may be following an eccentric orbit, as is the case for 
the majority of the Galactic globular cluster population, in which case the space velocity 
of Pyxis is unlikely to be perpendicular to its radius vector. There exists a set of proper 
motions for Pyxis (see Figure 4) that, when combined with V m , produce a space motion 
with a pole at (163, —22)° and leave the cluster bound to the Milky Way. 
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Fig. 4. — A plot of possible proper motions in equatorial coordinates for the Pyxis globular 
cluster. The interior of the circle defines the region of values for the transverse velocity 
that gives Pyxis a total space velocity less than the escape velocity from the Galaxy (~ 415 
km/sec). The shaded region defines the permissible proper motions that result in Pyxis 
having an orbital pole within 15° of the orbital pole of the LMC If Pyxis has been recently 
captured by the Milky Way, its proper motion should lie in this region. Given a proper 
motion in the shaded region and our measured radial velocity, the minimum eccentricity of 
the possible orbits will vary from ~ 0.70 — 1.0. The greyscale indicates the range of minimum 
eccentricities of the orbits given a proper motion in the various shaded areas. The scale is 
given by the bar below the plot. 
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Fig. 5. — The distribution of orbital eccentricities for globular clusters in the sample of 
Dinescu et al. (1999)| . The open histogram represents the eccentricities for the entire sample 
of 38 globular clusters with known orbits. The hatched histogram represents the distribution 
of eccentricities for the 10 objects in the sample that have R apo > 20 kpc, i.e., the outer halo 
objects. Given our large measured radial velocity for Pyxis, the only allowed orbits it may 
follow that nearly share the orbital pole of the LMC have e ~0.8. Since the majority of the 
outer halo clusters in the sample of pinescu et al. (1999)| have 0.6 < e < 0.8, it is plausible 



to assume that Pyxis may also follow such an eccentric orbit. 
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Table 1. Observed Pyxis Giant Stars 



ID 


^2000.0 


"2000.0 


# of Obs. 




R a 


(B - R) a 




(h:m:s) 


(d:m:s) 




(km/sec) 






Pyxis A 


09:08:01.7 


-37:13:48 


4 


32.7 


17.08 


2.01 


Pyxis B 


09:07:58.6 


-37:14:42 


2 


38.4 


17.75 


1.80 


Pyxis C 


09:07:56.2 


-37:13:16 


4 


36.6 


18.28 


1.78 


Pyxis D 


09:07:54.5 


-37:14:01 


5 


26.1 


18.33 


1.75 


Pyxis E 


09:07:57.8 


-37:13:17 


4 


37.9 


18.08 


1.78 


Pyxis F 


09:07:55.7 


-37:14:02 


3 


33.9 


18.09 


1.71 



Photometry from 
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